Load MG1655 genbank and generate a Karyotype
In [1]:
from pyrcos import utils
from Bio import SeqIO
In [2]:
record = SeqIO.read("../data/MG1655.gb", "gb")
karyotype = utils.seq_records_to_karyotype([record])
In [3]:
print(karyotype)
Load tracks from record
In [4]:
track = next(utils.seq_record_to_tiles([record]))
Create circos objects
In [6]:
from pyrcos.objects import Circos, Tile, Background
In [7]:
bg = Background(color="lgrey", y0=0, y1=1)
tile = Tile(file=track, r0=0.65, r1=0.95, backgrounds=bg, color="green", layers=5,
orientation="center", margin="0u", padding=50, thickness=50)
circos = Circos(karyotype, plots=tile, circos_path="/Users/joao/Documents/circos-0.67-7/", width=2500)
circos.attributes["max_points_per_track*"] = 5000
In [8]:
circos
Out[8]:
In [ ]: